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Abstract 

We consider the spin-i XY chain in a transverse field with regularly varying exchange interactions 
and on-site fields. In two limiting cases of the isotropic (XX) and extremely anisotropic (Ising) exchange 
interaction the thermodynamic quantities are calculated rigorously with the help of continued fractions. 
We discuss peculiarities of the low-temperature magnetic properties and a possibility of the spin-Peierls 
instability. 

PACS numbers: 75.10.-b 

Keywords: XX chain, Ising chain in a transverse field, regularly alternating spin chain 

It is generally known that there are a number of one-dimensional spin models for which in the translationally 
invariant (uniform) case an exact calculation of the thermodynamic quantities is possible. To this class belong 
the spin-i XY models proposed by Lieb, Schultz and Mattis [|J. At present there are relatively large number 
of magnetic and ferroelectric compounds which may be described by these models. On the other hand, due 
to the progress in material sciences a study of the properties of quantum spin chains with regularly varying 
parameters (e.g., exchange interactions or on-site fields) attracts much interest during last years. As an example 
one may mention the layered materials with the well-pronounced one-dimensional character in their properties 
known as superlattices, which can be viewed as chains with regularly varying exchange interactions (see, e.g., 
). In what follows we show how the thermodynamic properties of the regularly alternating one-dimensional 
spin-i XY model in a transverse field can be examined rigorously in two limiting cases of the isotropic exchange 
interaction (transverse XX chain) and of the extremely anisotropic exchange interaction (transverse Ising chain). 
We discuss briefly the effects of regular alternation on the magnetic properties of the chains at zero temperature. 
We also comment on a possibility of the spin-Peierls instability in the considered chains. 
We consider a ring of N — > oo spins \ described by the Hamiltonian 

h = J2 n»< + 2 CX< +1 + (i) 

n n 

concentrating on two limiting cases 1) of the isotropic exchange interaction 1% = 1% = I n (transverse XX chain) 
and 2) of the extremely anisotropic exchange interaction = J„, 1% = (transverse Ising chain). We assume 
that the exchange interactions /„ and the on-site (transverse) fields £l n vary regularly along the chain with 
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period p (N is a multiple of p), i.e., ■ ■ • ij,fipiifJii2^2 ■ ■ ■ I p Q P ■ ■ ■ ■ To calculate the thermodynamic 

quantities of the spin model we map it onto a system of spinless ( Jordan- Wigner) fermions with the Hamiltonian 



H = ^O n | <-;,<■„ 

n 

n ^ ' 
The Hamiltonian (2) can be brought into the diagonal form H = ^ fc A*, (77^% — \~) where is the elementary 
excitation energy and r/^ , r\k are Fermi operators. Obviously, the diagonalization of the Hamiltonian for p > 1 
becomes not trivial. 

To find the thermodynamic quantities of the spin model it is sufficient to know either p(E) = 5(E— A&) 

or R{E 2 ) — -k J2k 3(E 2 ~ ^1) since the Helmholtz free energy per site can be written as follows 

f3E\ 



f = -jjj dEp(E)\n ^2cosh^^ = -|jf dEER(E 2 ) In ^2 cosh ^ 



(3) 



Consider first the transverse XX chain. We introduce the one-fermion Green functions G nm — G nm (E) 
which satisfy the set of equations 

{E Q n *)G nm I n —\G n —\ ra InGn+l^rn — &nrn (4) 

and yield the required density of states p(E) = J2n I m G>m> GJ m = G nm (E±ie), e — > +0. One immediately 
notes that the continued fraction representation for the diagonal Green functions which follows from (4) 

E — Q„ — A n — An 
A" = A+ - — 



j2 i — n j2 

E - O n _! - K _o"J o 2 _ — E - tt n+ i - K _r " 



is extremely useful in the case of regularly varying parameters /„, f2 n since A~, A+ in (5) become periodic and 
can be evaluated exactly by solving quadratic equations. 

We cannot proceed in the described manner for the transverse Ising chain because of terms c+c^ +1 , c n c n +i 
in (2). However, from |lj we know that to find the (real) coefficients gu n = \ i^kn + ^kn), h>k n = \ {^kn — ^kn) 
of the canonical transformation r\k = Y) n {gknC n + hk n ct,) which diagonalizes the Hamiltonian (2) one should 
solve the set of equations 

^n-l4-l$fc,n-l + (^n + 7 n-l - A fc) ®kn + tt n /„$ fe ,„+i = (6) 

(and a similar set of equations for ^>kn)- Eq. (6) coincides with a set of equations describing displacements of 
particles in a nonuniform harmonic chain with nearest neighbour interactions and A& plays a role of oscillator 
frequency. We can find the distribution of oscillator frequencies squares R(E 2 ) from the relation R(E 2 ) = 
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^2 n Imr^„, r^ m = T nm (E 2 ± ie), e — > +0 where the Green functions T„ m = T nm (E 2 ) satisfy the set of 
equations 

From (7) the following representation for T nn may be derived 



r, m - 2 2 2 + , (8) 

^ "71 J n-1 *-»»» 

O 2 A 2 
_ "n-l J n -l 



o 2 r 2 

P2 _ O 2 T 2 "ti-2-'ti-2 

^ "71-1 1 n-2 £ 2 -SP „-/ 2 I 



fl 2 I 2 

71 77 



o 2 r 2 

^ "ti+1 J 7i B 2 -J2 2 , „-/ 2 , 



Again for a periodic sequence of parameters /„, f2„ the continued fractions A~, A+ in (8) become periodic and 
can be evaluated exactly by solving quadratic equations. 

To illustrate how the described approach works we consider the chains of period 1 (uniform chains). Then 
for the transverse XX chain one finds 



A- = A+ = ^[E-n±J(E-nf-u 2 ), 



-., a(E- ny < 4i 2 



p(E) = { * ^~(E-n-2i)(E-n+2i)' ' 
0, otherwise. 



For the transverse Ising chain the corresponding results read 



A" = A+ = 1 [E 2 -Q 2 -I 2 ± \]{E 2 -n 2 -I 2 ) 2 -4ft 2 / 2 

, -—j=, = f, if(E 2 -n 2 -i 2 ) 2 <4f> 2 / 2 , 

R{E 2 ) = \ V-( £2 -( n - 7 ) 2 )( s2 -( n + 7 ) 2 ) (10) 
0, otherwise. 

Formulas (3), (9), (10) reproduce the long known results, however, the described approach (formulas (3), (5), (8)) 
yields new results for the chains of longer periods. These calculations are more lengthy (although transparent) 
and therefore are omitted here (see ||, [|] ) . 

Let us discuss the effects caused by regular nonuniformity which can be studied rigorously within the 
framework of the elaborated approach. We start from the transverse XX chain. Regular alternation of the 
Hamiltonian parameters leads to a splitting of the fermion band (see (2)) into several subbands. Since the 
(transverse) magnetization per site is m = — ^ J!^, dEp(E) tanh 4p one immediately observes the plateaus in 
the dependence m vs. fl at zero temperature (3 — > oo (we assume Q n = f2 + Af2„) which are accompanied 
by the square-root singularities in the zero-temperature dependence (transverse) susceptibility \ vs - ^- The 
number of plateaus does not exceed the period of alternation, the characteristic fields at which the plateaus 



3 



start and end up (and the singularities of \ occur) can be easily calculated from the density of states p(E). 
Consider further the case of modulated exchange interactions \Ii^\ = ± 8) (0 < 5 < 1 is the dimerization 
parameter) and constant field f2. One may easily convinced himself that the ground-state energy per site 
eo = — \ dEp(E)\E\ for |fi| < 2\I\ decreases sufficiently rapid with increasing of 5 to provide a spin-Peierls 
dimerization. Moreover, the ground-state energy of the chains having longer periods of exchange interaction 
modulation obtained with the help of continued fractions allows one to examine the stability of the transverse 
XX chain with respect to more complicated lattice distortions. 

Let us pass to the transverse Ising chain. Considering the case of modulated exchange interactions (e.g., 
of period 2) and constant field and using the strong-coupling approach arguments (i.e., assuming the smallest 
interaction to be equal to 0) H one immediately notes that no magnetization plateaus should be expected. 
However, using the continued fraction approach results one may show that a logarithmic singularity of x may 
occur if Q = ±1 (p — 1), = H1I2 (p = 2) etc. As a result the transverse Ising chain may exhibit 

'plateaus' in the dependence m vs. O for sufficiently large deviations of fields Af2 n . Moreover, in contrast to 
XX chain, the number of singularities of x (and the number of 'plateaus') for given p depends on the values 
of Hamiltonian parameters. The direct calculation of the ground-state energy of the dimcrized Ising chain in a 
constant transverse field shows that no spin-Peierls instability should be anticipated that is in agreement with 
general arguments about the Peierls mechanism. Note, however, that the transverse Ising chain is unitarily 
equivalent to the anisotropic XY chain without field ||. The latter model for a special choice of parameters 
may become the isotropic XY model for which, as we know from the rigorous calculation, the spin-Peierls 
instability does occur. 

To summarize, we have discussed an application of continued fractions for calculation of the thermodynamic 
quantities of the regularly alternating spin-i transverse XX and Ising chains. We have briefly discussed the 
magnetic properties of such chains as well as their stability with respect to a lattice distortion due to the 
spin-Peierls mechanism. 
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